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Abstract 

We examine the radiative corrections to the mass of the hghtest Higgs 
boson in the minimal supersymmetric extension of the standard modeL 
We use the renormahzation-group improved effective potential which in- 
cludes the next-to-leading-order contributions. We find that, contrary 
to the result of Espinosa and Quiros, the higher-order corrections to the 
lightest Higgs boson mass are non- negligible, adding 3 — 11 GeV (3 — 9 
GeV) to the result in the leading logarithm approximation for the range 
of top quark mass lOOGeV < < 200GeV and for the supersymmetric 
breaking scale MgusY = ITeV (Msusy = lOTeV). Also we find that our 
result is stable under the change of the renormalization parameter t. 
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Although the standard model (SM) is highly successful and in excellent agree- 
ment with all the measurements at the present energies, it is widely believed that 
SM is not the final theory for the world of elementary particles. The minimal super- 
symmetric extention of the standard model (MSSM) is one of the most promising 
candidates beyond the SM. The MSSM possesses in its physical spectrum three 
neutral and two charged Higgs bosons, and there exists a tree-level relation which 
implies that at least one neutral Higgs boson is lighter than the mass [Mz)- 
Radiative corrections to the masses of these Higgs bosons have been calculated 
by several groups |l|] - who found that they are quite significant, depending 
strongly on the top quark mass and the scale of supersymmetry breaking (Msusy) 
or the squark masses. All the above works except ref. considered the one-loop 
radiative corrections to Higgs boson masses. 

Now that the one-loop corrections have been found to be significantly large, it is 
quite natural to ask next how large the higher-order corrections would be. Indeed 
Espinosa and Quiros |]^ have analyzed the "two-loop" radiative corrections to the 
mass of the lightest Higgs boson in the minimal and non-minimal (including a gauge 
singlet) supersymmetric standard model. They used the effective potential (EP) 
in the leading logarithm approximation and examined the evolution of the Higgs 
quartic coupling A by the renormalization group (RG) techniques with the one- and 
two-loop /? functions. They found that the "two-loop" correction is negative and 
stays within a few percent even in cases where the one-loop correction is larger than 
the tree-level mass. 



Recently there appeared interesting papers [15| - |17| which discussed about 



the improvement of EP by using the renormalization group equation (RGE) . It 
was shown there that to improve EP which satisfies the RGE with up to the two- 
loop /5 functions and anomalous dimension 7, one should include the one-loop- 
level potential with the running parameters into the solution. In this respect, the 
work of Espinosa and Quiros seems unsatisfactory: they used the EP in the 
leading logarithm approximation, which is the tree-level potential with the running 
parameters, and they made use of the two-loop /3 functions only to determine the 
evolution of these parameters. In this paper we reanalyze the mass of the lightest 
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Higgs boson (m^) in the MSSM using the EP improved by RGE up to the next-to- 
leading order. We find that new terms which were not considered by Espinosa and 
Quiros give non-neghgible contributions to the m^. We also find that the predicted 
values of are stable under the change of the renormalization parameter t when we 
use the RGE-improved EP which includes the next-to-leading-order contributions. 

Two basic assumptions were made in their analysis of the lightest Higgs boson 
mass 0: (a) all supersymmetric (SUSY) partners of the SM particles have masses 
of the order of the supersymmetry breaking scale Msusy; (b) one linear combination 
H of the two Higgs boson doublets, = {H^, H{f and H2 = {H^, H^f, 

H = H^cosp + iT2H;smp (1) 

is light, while the other linear combination, which is orthogonal to the former one, 
is as heavy as the SUSY partners. Under these assumptions, it is clear that the 
effective theory below the scale MgusY is the usual SM with one light Higgs doublet 
H. Throughout the following analyses, we will make the same assumptions (a) and 
(b). The tree-level Higgs potential below Msusy is then given by 

V^tree = -m2|i^|2 + lA|i^|^ (2) 


where 

lx = \i9! + 9'2)cos'2P, (3) 

and gi and g2 are the gauge coupling constants of U{1)y and SU{2)l, respectively. 

When the neutral component of H acquires a vacuum expectation value < > = 
vj -\/2, the above tree- level potential (H) gives the physical Higgs = (Rei7° — 1>)/-\/2 
(which corresponds to the lightest Higgs boson under the assumptions (a) and (b)) 
a squared mass 

m\ = Ixv'. (4) 

Also at the symmetry breaking, the top quark and the Z^, gauge bosons gain 
masses which are given by 

mt = htv/V2, Ml = \{gl + gl)v\ = -^glv\ (5) 
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where ht is the Yukawa couphng of to the top quark. The tree-level relation 



m 
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M|cos22/5 follows from Eqs. (D - (D- 
The (RG15-unimproved) EP of the SM up to the one-loop level is given by 

V, = 1^(0) + ^(1), (6) 
^(0) = -^mV + ^A^^, (7) 

where (pc is the classical field corresponding to the physical Higgs boson 0, and all 
the Yukawa couphngs of to quarks and leptons except the top quark are neglected. 



The calculation is performed in the Landau gauge and in the MS scheme to obtain 
the one-loop result V(i), and /i is the renormalization scale. The ellipses in Eq. (H) 
represent contributions of the gauge bosons and the would-be-Goldstone bosons. 
Throughout this paper we use the Landau gauge which is the most convenient for 
our purpose 



THl and the MS scheme. 



Now we improve the EP by using the RGE. It was recently emphasized by the 



authors of Refs. |]TB| - [|T^ that in the MS scheme the EP V{(pc) fails to satisfy the 
usual (homogeneous) RGE unless V^(0), a contribution to the "vacuum energy", is 
suitably dealt with. When we use the RGE-improved EP in the leading order and 
obtain m^, the consideration of V{0) term is unnecessary. However, as we shall see 
below, if we improve the EP by RGE up to the next-to-leading order, V{0) becomes 
relevant to us and we must take its presence into account. Thus with an appropriate 
(y9c-independent term being added, the new EP V{(pc) satisfies the following RGE of 
the usual form: 

{V-j^^,—)V{y^„Xi,fi) = (9) 

with 

where Xj = A, ht, m^, g^, g2, gi and g^ is the gauge coupling constant of SU {S)c- The 
solution is easily found by the method of characteristics and we get 

V{if,,X,,fi) = VMt),X,{t),^{t)), (11) 
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m 


= exp | — 








= /ie*. 



where Xi{t) are running couplings and running mass which are determined by the 
equations 

= PxA^M, X„Xj = X,ht,m',gs, 92,91 (12) 
with the boundary conditions Xj(0) = Xi , and 

(13) 



Then using the result of tree- and one- loop-level EP of Eqs.(^ and (||), we obtain 
the RGE-improved V as follows: 

V = Q{X,{t),fi{t)) 

+V^o){Vc{t),Xi{t)) 

+V(i)(v9e(t),X,(t),/x(t)) + ---, (14) 

where Q is the y^c-independent term which is added for V to satisfy a RGE of the 
usual homogeneous form, and the ellipses represent the higher-loop contributions. 

For later convenience, let us expand the RGE coefficient functions j3xi and 
by the number of loops as follows: 

(^x, = ^f^x] + ^'^(^x] + • • • , Xi = \ht, ^3, 92, 91 

14> = H^+h'lf + ---. (15) 

where we have introduced the Plank's constant h so that the power of h counts 
the number of loops and l3x] and 7^"^ are the n-loop contribution to jSx^ and 7^, 
respectively. Similarly V has the loop expantion 

V = ^ + V(o)(v5c(t)) + hV^i){Vc{t)) + ■ ■ ■ , (16) 

and we have denoted i7(Xj(t), /i(t)) and V(„)(v9c(t), Xj()f:), /i(t)) as VL and V{n){ipcit)), 
for short, respectively. 

Inserting Eq.(0) into Eq.(^ and picking the terms up to of the order Ti, we find 

.„,a{,«^-,<V^.JL^.o, ,17) 
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where the last term in the parentheses arises from ^dV{i)/d^ when we use the 
expression of V(i) in Eq.(||) and neglect the contributions of the gauge bosons and 
the would-be- Goldstone bosons to V{iy If we further set = in the above 
equation, we obtain 

p^] + ;.^^^i#^ = o, (18) 

from which Vt can be determined to the leading order. For the later discussion, 
however, we do not need the specific form of Vt , and we will see below that the 
knowledge of Eq.([T7|) is sufficient for our purpose. 

We will now analyze the mass of the lightest Higgs boson using the RGE- 
improved V . At first, let us consider the boundary conditions for coupling constants. 
Under our basic assumptions explained before, the relation between the quartic cou- 
pling constant A and the gauge coupling constants gi and g2 given in Eq.(^) should 
be satisfied at the scale Msusy, that is, 

^A(Msusy) = \{gl{MsvsY) + (?2'(^susy))cos'2/5. (19) 

So we choose the renormalization scale /i to be MgusY and take the parameter t as 
t = In (y9c/MgusY)- Then we find that /i(t) = and the RGE-improved V is given 
by 

V = n + r(o)(</^e(t)) + mi)iVc{t)) + 0{n^), (20) 

with 



F(i)(^,(t)) = -^^ht{t)^t{t) 



2 2 



(21) 



where in V(i){(pc{t)) we only include the top-quark contribution to the one-loop 
EP, because, with its very large Yukawa coupling, the contribution of top-quark is 
dominant over those from the gauge bosons and the would-be Goldstone bosons. We 
will obtain the lightest Higgs boson mass by evaluating ff^V /dipcitY at Pcitv) = v = 
{V2Gf)-2 = 246GeV under the minimum condition dV /d(fc{t) = at (fc{tv) = v. 
The value is determined by the equation (pdtv) = v, which, with help of Eq.(|T3[), 
is transformed into 

t,=\nj^+ f\^{t')dt'. (22) 



SUSY 
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It is noted that we differentiate V not by (fc but by the renormahzed field ipc{t) and 
also that we evaluate the differetials at the point of ipc{tv) = v and not at = f . 
Since f^cif)^ and are functions of t, we find 



/i(t) 



dip 



d(pc{t) ipcit) 



+ o{h). 



(23) 



Thus we obtain 



d(pc{t) 
dipcit) 



(24) 



-n- 



.2"-* 



In 



hmeit) 



It is now straightforward to evaluate dV/d(fc{t) from Eq.(|^). Using the relation 
(p!?!), we eliminate T>Q term and find 



dV 



+0(/i' 



167r- 



2 



(25) 



Further differentiation of dV/dipdt) by fcit) gives 



hmeit) 
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./.^(t)^^(t) 



In- 



+C»(^'). 



(26) 
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Using the minimum condition dV/d(fc{t) = at (fc{tv) = v, we finally obtain for 
the lightest Higgs mass in the next-to-leading logarithm approximation 



m 



</)(2-loop) 



(Pc{tv)=V 



(27) 



Stt- 



;ht(ty) 



In 



} 



The first term in Eq.(^) gives 



rrij 



;i-ioop) 



(1/3) x{t,y, 



(28) 



which is the result given by Refs. |^ except that the authors of Ref. have 
evaluated the running coupling X(t) at = ln(m0/MsusY)- As far as the leading 
logarithm approximation is concerned, the terms of order h in Eq. (|27|) are neglected 
as the higher-order effects. An arbitrariness coming from different choice of the 
parameter t also falls in the higher-order corrections although it has no small effect 
on the predicted values numerically. Espinosa and Quiros have employed the "one- 
loop" formula = {l/3)X{t)v^, which is correct only in the leading logarithm 
approximation, and computed \(t) at t^^ = ln(m(^/MsusY) using RGE with up to 
two- loop (3 functions 0. However, if we make use of the two- loop RGE coefficient 
functions for the running parameters and evaluate the lightest Higgs mass m^, we 
should take into account the order-Zi-terms in Eq. (pTj) which also collect the next- 
to-leading logarithmic contributions. 

The one- and two-loop (3 functions and anomalous dimension 7^ for the SM, 
which we will use in this analysis, read as follows |T^. We define the constant A as 
A = 167r2. 

For the Higgs quartic coupling A: 



+ 12A/i^ 



3Qht - SXiSgl + gl) + '^^gt + \glgl + \gl 
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27 57 
+180/1,6 - ht{192gl + 16gf) + hl{--gt + 63^2^^^ _ _^4) ^29) 

915 289 4 2 559 ^ . 379 « 
+ -^^2 - -^9291 - -^9291 - -^91- 

For the top-quark Yukawa coupling ht- 

= \hl-h,{Sgl + \g,+ ^lg% 



A'i^t} - h,{-12ht-2Xh',+h',{36gl+'-^g', + ^gl) (30) 
+ 1a^ - 108,3^ + 9glgl + '-^g^gl - |,,^ - ^glgl + ^,^}. 
For the gauge couphngs g^, g2, and gi. 

APS,^ = -7<?i 

^/^i^^ = -'{91 (31) 
APi'^ - 'i9h 

^'P?^ - gli-'^hl + ^^gl + lgl+'^gl). 
For the mass parameter m^: 

= m^(2A + 6/^?-^^i-^^D, 
^^/^S = m^{-^A^-12A/.,^ + 4A(3^,^ + ^?)-f/^^ (32) 
+hl{AQgl + ^gl + - ^(72' + {9I9I + ^^^}- 

For the anomalous dimension 7^: 
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l2^,(2) _ ^2 27 4 ,2/on^2 I 45 2 , 85 2x 



= ^y-^ht + ht{20gi + -g', + --gt) (33) 



6 4 ' 8 24 

271 4 , 9 2 2 431 4 



92 + T^^2^i + TT^^i- 



32 16"^"' 96 

When we substitute the expressions of the RGE coefficient functions l3^^\ (3^2, 
and 7^^'' into Eg . (1271) , we find that many terms cancel out and we obtain a rather 
simple expression 

2 _ ■'■ \ 2 

'^(/.(2-loop) — 2^"^ 

+^^"-[t~2{l^" + 2A/^? - ^^(3^2^ + gl) (34) 

where it is understood that all the running parameters are the ones evaluated at 
t = ty. It is interesting to note, in particular, that a term in Eg . (|27|) , which comes 
from the one-loop EP of Eg.(|^), cancels with another term in (3^^ of Eg.(p9|) 
and, in conseguence, there appears only one hf term of the form hfln{h'^C,'^ /2) in 
the order-^ contributions in Eg. (|3^ . This makes the higher-order corrections to 
m|(2-ioop) to be a milder one. From the expression of Eg.(0) we expect that the 
next-to-leading-order corrections give a positive contribution to m^, which will be 
shown numerically to be true below. 

To evaluate m(^(2-ioop) of Eg. (|3^ ) numerically, we choose the initial conditions for 
the gauge couplings ctj = gf/4:TT (i=3,2,l) at the scale Mz = 91.2GeV to be 

as{Mz) = 0.115, a2iMz) = 0.0336, ai(M^) = 0.0102 (35) 

which are consistent with present experimental constraints [|l^ - |2ll], and define the 
Yukawa coupling of the top guark at the scale of its mass rrit as 

ht{mt) = V2mt/v with v = 246GeV. (36) 

The fact that the two-loop (3^^ , and Z?^^-' (i=3,2,l) are functions of the couplings 
A, ht, and gi casts Eg.([T^) into a very complicated system of coupled differential 
eguations. For given values of mt and Msusy, we first solve the system (|12D and 
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( P^D with jSxi and 7^ given in Eqs.(pUD - ( P5| ) together with the initial conditions 
Eqs.(^5|) and (|36|) , and we obtain the appropriate and X(t^) so that A satisfies 
the boundary condition, Eg . (|19D , when it evolves from to t = 0. At the same 
time when we find the appropriate and A(tt,), we gain all the information on the 
parameters which appear in the r.h.s. of Eg . (|3^) . This is how we calculate m(^(2_ioop) 
for given values of rrit and Msusy- 

In Fig.l we plot m^(2_ioop) as a function of mj for MgusY = 1 TeV, cos^ 2/3 = 

1 and cos^ 2/5 = along with m^(i_ioop) in the leading logarithm approximation. 
Fig. 2 shows the case for Msusy = 10 TeV. Since it is suggested that rrit is not too 
excessively large in MSSM ||2^, we have studied m0(2-ioop) for rrit from lOOGeV 
to 200GeV. The curves for cos^ 2/3 = 1 (cos^ 2/3 = 0) can be considered as upper 
(lower) bounds for the lightest Higgs boson mass in the MSSM. From Fig.l and 

2 we observe that the next-to-leading-order effects are non- negligible. They add 

3 — 11 GeV (3 — 9 GeV) to the result in the leading logarithm approximation for 
the range of top guark mass lOOGeV < rrit < 200GeV and for Msusy = ITeV 
(Msusy = lOTeV). These rather large corrections come from the order-^-terms of 
Eg. (|3^) , since without those terms we could indeed recover the result of Ref. 0, 
namely, the higher-order corrections being negative and negligible for the considered 
range of parameters. Contrary to the conclusion of Espinosa and Quiros, our result 
shows that the higher-order corrections turn out to be positive and non-negligible 
when top guark is very heavy. 

It is to be noted that our result in the leading logarithm approximation differs 
numerically from those of Refs. 01 because we evaluated the running coupling 
\{t) at tv given by Eg. (p2|) , instead of at tm^ = ln(m<^/MsusY)- In other words, 
we have made a different choice of t from the ones made in the above references. 
In the case of cos^2/3 = 1, for example, our predicted values for m^(i_ioop) are 
smaller than those calculated in Refs. [|I| by — 6GeV (0 — lOGeV) for the 
range lOOGeV < rrit < 200GeV and for Msusy = ITeV (Msusy = lOTeV). In the 
leading logarithm approximation, the predicted values for the lightest Higgs boson 
mass are rather sensitive to the choice of the renormalization parameter t. However, 
in the next-to-leading logarithm approximation, it is not possible to change the 
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definition of t witliout modifying the 0{h) terms of EP. In consequence, the result 
is stable under the change of t. This is a well-known issue which often arises in the 
renormalization group approach, and will be discussed in more detail in the second 
comment below. 

Fig.l shows that the next-to-leading-order corrections are large for cos^2/3 = 
especially when Msusy = ITeV. In the case of cos^ 2/? = 0, the boundary condition 
for X{t) at Msusy is A(t = 0) = 0. When Msusy = ITeV, the "evolution time" is not 
long enough for X{t) to grow from the initial value A = 0, and thus a term of the form 
—6hfln{hfC,'^/2) and gauge-coupling-constant terms in the order-^ contributions in 
Eq.(|3^) give relatively large corrections compared with the leading (l/3)A(tt,)f^. 

A few comments are in order. The first comment concerns the definition of mass. 
The "mass" we have calculated is not the on-shell mass. It might be necessary to 
consider the correction coming from the wave function renormalization in order to 
make a realistic prediction for the Higgs mass. However, this effect is expected to 
be small § 

Secondly, we have taken the parameter t as t = ln((y9c/MsusY) to derive Eq.(p7|). 
But the physics should not depend on the choice of t. For example, as was stressed 
in Refs. [|l^ , the "natural choice" of t may be given by the equation 

2^^\t) = 2^^'e'' = hl{t)vl{t). (37) 

Then the RGE-improved EP which we deal with will be 

V = fi(x,(t),Mt))-^m'Wv^^W + ^A(t)v.^(t) 

+^{-[^2htit)Vc{t)]. (38) 

Since the expression of the order-^ term of the above equation is different from 
the one we analysed before, we may think at first sight that we would obtain a 
different result for m(^(2-ioop)- In fact we follow the same procedure as before, i.e., 
first differentiate V by fc{t), use the relation (|17D and eliminate term, evaluate 
d^V /df^city at ipc{tv) = V under the minimum condition dV/d(pc{t) = at ipdtv) = 
V, and we obtain the following expression: 

^(2-ioop) = ^HQv'^ 
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(39) 



87r2 



A logarithmic term of the form —{3/8'K^)hv'^hfln{hfC,'^/2) which appeared in Eg . (1271) 
is missing. However, we should note that the definition of t has been altered. Re- 
membering that we have chosen the renormalization scale /i to be Msusy, we must 
now evaluate the curvature of EP at the value of ty, which is determined by 



In- 



+ In 



SUSY 



ht{tv) 



(40) 



Expanding \{tv) around t^, we find 



A(t.) 



^7! ti 



+ 



(41) 



Then substituting the above expression for A(t^) into Eg. (p9D , we will obtain essen- 
tially the same result as Eg. ( p7D except for the h"^- and other non-dominant terms. In 
fact, we have calculated the mass m<^(2-ioop) with the above choice of t, Eg.(0), and 
found the difference between the two numerical results for m(^(2_ioop) and m0(2_ioop) 
being less than IGeV for < 200GeV and MgusY < lOTeV. 

In conclusion we have examined the mass of the lightest Higgs boson in the 
MSSM beyond the leading logarithm approximation. We have made use of the EP 
improved by RGE up to the next-to- leading order. We have found that, contrary to 
the result of Espinosa and Quiros, the next-to-leading-order corrections to the Higgs 
mass are non-negligible, adding 3 — llGeV (3 — 9GeV) to the values predicted by the 
RGE approach in the leading logarithm approximation for the range lOOGeV < rrit < 
200GeV and for MgusY = ITeV (MgusY = lOTeV). We also found that the predicted 
values of are stable under the change of the renormalization parameter t when we 
use the RGE-improved EP which includes the next-to-leading-order contributions. 
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Figure caption 

Fig.l 

Values of the Higgs boson mass as a function of rrit, for Mgusv = ITeV and cos^ 2/3 = 
1 (upper two lines) and cos^2/9 = (lower two lines). The sohd and dash-dotted 
lines denote the next-to-leading-order and the leading-order results, respectively. 

Fig.2 

The same as in Fig. 1, but considering the case for Msusy = lOTeV. 
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